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The leading-log order QCD hard scattering Monte-Carlo models of HERWIG, ISAJET, and 
PYTHIA are used to study the sources of &-quarks at the Tevatron. The reactions responsible 
for producing b and b quarks are separated into three categories; flavor creation, flavor excitation, 
and parton-shower/fragmentation. Flavor creation corresponds to the production of a bb pair by 
gluon fusion or by annihilation of light quarks, while flavor excitation corresponds to a fe or b quark 
being knocked out of the initial-state by a gluon or a light quark or antiquark. The third source 
occurs when a bb pair is produced within a parton-shower or during the fragmentation process of 
a gluon or a light quark or antiquark (includes gluon splitting). The QCD Monte-Carlo models 
indicate that all three sources of 6-quarks are important at the Tevatron and when combined they 
qualitatively describe the inclusive cross-section data. Correlations between the b and b quark are 
very different for the three sources and can be used to isolate the individual contributions. 



I. INTRODUCTION 

It is important to have good leading order (or leading-log order) estimates of hadron-hadron collider observables. 
Of course, precise comparisons with data require beyond leading order calculations. If the leading order estimates are 
within a factor of two of the data, higher order calculations might be expected to improve the agreement. On the 
other hand, if the leading order estimates are off by more than about a factor of two of the data, one cannot expect 
higher order calculations to improve the situation. In this case, even if the higher order corrections were large enough 
to bring agreement, one could not trust a perturbative series in which the second term is greater than the first. If a 
leading order estimate is off by more than a factor of two, it usually means that one has overlooked some important 
physics. For this reason good leading-log order estimates are important. 
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FIG. 1: Illustration of a proton-antiproton collision in which the QCD hard 2-to-2 reaction results in the creation of a bb pair 
via the subprocess q + q ^ b + b or g -\- g ^ b -\-b. The event also contains possible initial and flnal-state radiation and particles 
that result from the break-up of the initial proton and antiproton (i.e., "beam-beam remnants"). 

In this analysis the leading- log order QCD hard scattering Monte-Carlo models of HERWIG |jl| , ISAJET |^ , and 
PYTHIA are used to study the sources of 6-quarks at the Tevatron. The reactions responsible of producing 
6-quarks are separated into three categories; flavor creation, flavor excitation, and shower/fragmentation. Flavor 
creation corresponds to the production of a bb pair by gluon fusion or by the annihilation of light quarks via the 
two 2-to-2 parton subprocesses, q + q ^ b + b, and g + g — > b + b, and is illustrated in FIG. The data from 
CDF and DO ^, ^, ^, Q for the integrated inclusive 6-quark cross section for \y\ < 1 at 1.8 TeV are compared with 
the QCD Monte-Carlo model predictions for flavor creation in FIG. 0, where y is the rapidity of the 5-quark. Here 
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FIG. 2: Data from CDF and DO jj, g, Q, [7| for the integrated inclusive 6-quark total cross section (pt > pT(min), \y\ < 1) 
for proton-antiproton collisions at 1.8 TeV compared with the QCD Monte-Carlo model predictions of HERWIG 5.9, PYTHIA 
6.115, and ISAJET 7.32 for the "flavor creation" subprocesses illustrated in FIG. ^ All three Monte-Carlo models were 
generated using the parton distribution functions CTEQ3L and pT(hard) > 3GeV/c. 



the parton distribution functions CTEQ3L have been used for all three Monte-Carlo models and, as is well know, 
the leading order predictions are roughly a factor of four below the data. The leading order estimates of the flavor 
creation contribution to 6-quark production at the Tevatron are so far below the data that higher order corrections 
(even though they may be important) cannot be the whole story. 




FIG. 3: Illustration of a proton-antiproton collision in which the QCD hard 2-to-2 reaction corresponds to the scattering of a 
6-quark out of the initial-state into the final-state by a gluon or a light quark or light antiquark via the subprocess g + b g + b, 
q + b q + b, or q + b q + b. These subprocesses together with the corresponding b terms are referred to as "flavor excitation". 
The event also contains possible initial and final-state radiation and particles that result from the break-up of the initial proton 
and antiproton (i.e., "beam-beam remnants"). 

An additional source of 5-quarks at the Tevatron comes from the scattering of a 6 or 6 quark out of the initial-state 
into the final-state by a gluon or by a light quark or antiquark via the subprocesses; g + b—>■g + b,g + b—^g + b, 
q + b ^ q + b, q + b ^ q + b, q + b q + b, and q + b ^ q + b. This is referred to as "flavor excitation" and is illustrated 
in FIG. ^. Flavor excitation is, of course, very sensitive to the number of b and b quarks within the proton (i.e., the 
structure functions). The b and b quarks are generated through the evolution of the structure functions. Even 
with no "intrinsic" 66 pairs within the proton, at high bb pairs are produced by gluons and populate the proton 
"sea" . The number of bb pairs within the proton is related, through the evolution, to the gluon distribution within 
the proton. None of the structure functions considered in this analysis include "intrinsic" 66 pairs within the proton. 
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The bb pair content within the proton is generated entirely through the evolution of the structure functions. 




FIG. 4: Illustration of a proton-antiproton collision in which a bb pair is created within a parton-shower or during the frag- 
mentation process of a gluon or a light quark or antiquark. Here the QCD hard 2-to-2 subprocess involves only gluons and 
light quarks and antiquarks (no heavy quarks in the 2-to-2 hard scattering subprocess) . The event also contains possible initial 
and final-state radiation and particles that result from the break-up of the initial proton and antiproton (i.e., "beam-beam 
remnants" ) . 

Another source of 6-quarks at the Tevatron comes from reactions which have a bb in the final-state but with only 
gluons and light quarks and light antiquarks participating in the 2-to-2 hard parton scattering subprocess (i.e., no 
heavy quarks in the 2-to-2 hard scattering subprocess). This is referred to as "shower/fragmentation" and is illustrated 
in FIG. 0. Here the subprocesses are all QCD 2-to-2 gluon, light quark, and light antiquark subprocesses. The 
"shower/fragmentation" contribution comes from bb pairs produced within parton-showers or during the fragmentation 
process. This category includes the "gluon splitting" subprocess, g + g ^ g + {g ^ bb), as modeled by the QCD 
leading-log Monte-Carlo models. 

Section II discusses methods of generating the three sources of 6-quarks; flavor creation, flavor excitation, and 
shower/fragmentation. In Section III the leading- log QCD Monte-Carlo model predictions are compared with data on 
the 6-quark inclusive cross section at the Tevatron. Correlations between the b and b quark are very different for the 
three sources and in Section IV we study ways to isolate the individual contributions using bb correlations. Section V 
is reserved for summary and conclusions. 

II. MONTE-CARLO GENERATION 

It is not easy to arrive at the QCD Monte-Carlo model predictions for 6-quark production at the Tevatron. Some 
of the plots presented here contain 8 million generated events. In addition, one must handle each of the Monte-Carlo 
models differently in order to get all three contributions; flavor creation, flavor excitation, and shower/fragmentation. 

A. Flavor Creation 

The flavor creation contribution to 6-quark production at the Tevatron is the easiest to generate. ISAJET, HERWIG, 
and PYTHIA all allow the user to select and generate separately the flavor creation contribution. For ISAJET if 
one runs all QCD and then selects the flavor creation terms by monitoring the 2-to-2 subprocess, one gets the same 
answer as running the flavor creation contribution separately. However, this is not true for HERWIG and PYTHIA. 
These two Monte-Carlo models only include the 6-quark mass when one runs the flavor creation terms separately 
{i.e., heavy quark production). For HERWIG and PYTHIA, when one runs all QCD the 6-quark mass is set to zero. 
ISAJET always includes the 6-quark mass in the flavor creation subprocesses. Setting the 6-quark mass equal to zero 
increases the fe-quark inclusive flavor creation cross section for pT> 5 GeV/c and \y\ < 1 by roughly 40% at 1.8 TeV. 

In summary, one can produce the flavor creation contribution in two ways with ISAJET, either by running them 
separately or by running all QCD. For HERWIG and PYTHIA one must generate the flavor creation terms separately 
by running the heavy quark option (IPROC = 1705 or MSEL ~ 5). For PYTHIA this produces only the flavor 
creation contribution. However, for HERWIG the heavy quark option produces both the flavor creation terms and 
the flavor excitation terms and one must monitor the 2-to-2 subprocesses to separate them. 
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FIG. 5: Data on the integrated inclusive &-quark total cross section (pr > PT(min), \y\ < 1) for proton-antiproton collisions at 
1.8 TeV compared with the QCD Monte-Carlo model predictions of PYTHIA 6.115 (CTEQ3L, pT(hard) > 3 GeV/c). The four 
curves correspond to the contribution from flavor creation (FIG. 0), flavor excitation (FIG.]^, shower/fragmentation (FIG. 
and the resulting total. 



B. Flavor Excitation 



For ISAJET one can run the flavor excitation contribution separately or one can run all QCD and then select out 
the flavor excitation contribution (one gets the same answer either way). For HERWIG running the heavy quark 
option gets both the flavor creation and the flavor excitation contributions (with the 5-quark mass included). For 
HERWIG, if one runs all QCD and then selects out the flavor excitation contribution one gets a different answer since 
here the b mass is set to zero. For PYTHIA the only way to get the flavor excitation contribution is to run all QCD 
and then select out the flavor excitation contribution, but here one only gets the massless approximation. Actually, 
in leading order it is not clear how to handle the 6-quark mass in the flavor excitation subprocesses since one is going 
from an off-shell b (or b) quark in the initial state to an on-shell massive b (or b) quark in the final state. This limits 
the accuracy of the leading order estimate of this contribution. 



C. Shower /Fragmentation 



For all three Monte-Carlo models the only way to get the shower/fragmentation contribution to 6-quark production 
is to run all QCD 2-to-2 subprocesses and then select out the shower/fragmentation contribution by monitoring 
the subprocesses, which requires producing large numbers of events. The shower/fragmentation contribution differs 
greatly among the Monte-Carlo models for two reasons. The first is due to different fragmentation schemes. ISAJET 
uses independent fragmentation, while HERWIG and PYTHIA do not. The second difference arises from the way the 
QCD Monte-Carlo produce parton showers. ISAJET uses a leading-log picture in which the partons within the shower 
are ordered according to their invariant mass. Kinematics requires that the invariant mass of daughter partons are 
less than the invariant mass of the parent. HERWIG and PYTHIA modify the leading- log picture to include "color 
coherence effects" which leads to "angle ordering" within the parton shower. 

The flavor excitation and shower/fragmentation 2-to-2 subprocesses diverge as the parton transverse momentum, 
PT(hard), goes to zero. To get a finite result one must regulate the divergences or specify a minimum value for 
PT(hard). For HERWIG and ISAJET we take ]3T(hard) > 3GeV/c. PYTHIA regulates the perturbative 2-to-2 
parton scattering differential cross sections, which allow one to take pT(hard)> 0. For PYTHIA 6-quark predictions 
with pT(fe-quark) > 5GeV/c are independent of the whether one takes pT(hard) > 3GeV/c or pT(hard)> 0. 
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FIG. 6: Data on the integrated inclusive 6-quark total cross section (pr > pT(min), \y\ < 1) for proton-antiproton collisions at 
1.8 TeV compared with the QCD Monte-Carlo model predictions of ISAJET 7.32 (CTEQ3L, pT(hard) > 3GeV/c). The four 
curves correspond to the contribution from flavor creation (FIG. 0), flavor excitation (FIG. shower/fragmentation (FIG. 
and the resulting total. 



III. THE INCLUSIVE CROSS-SECTION 

The data from CDF and DO ||, ||, 0] on the integrated inclusive &-quark total cross-section are compared with the 
QCD Monte-Carlo model predictions in FIGS. ^ with the CTEQ3L structure functions and pj{ha,Td) > 3GeV/c. 
The four curves in each of the plots correspond to the contribution to 6-quark production from flavor creation, 
flavor excitation, shower/fragmentation, and the resulting overall total. After adding the contributions from all three 
sources PHYTHIA (CTEQ3L) agrees fairly well with the data. ISAJET (CTEQ3L) is a bit below the data because 
of a smaller shower/fragmentation contribution. HERWIG (CTEQ3L) is also slightly below the data due to a smaller 
flavor excitation component. One does not expect precise agreement from these leading-log estimates. However, the 
qualitative agreement indicates that probably nothing unusual is happening in 6-quark production at the Tevatron. 

FIG. 1^ and FIG. || show the differential inclusive 5-quark px cross section and the integrated total cross section, 
respectively, from PYTHIA (CTEQ4L), with pT{ha,rd)> 0. PYTHIA regulates the divergences in the perturbative 
2-to-2 parton scattering differential cross sections and these plots show a smooth behavior down to pT(&-quark) equal 
to zero. 

FIG. M shows the predictions of HERWIG (DOl.l), HERWIG (CTEQ3L), ISAJET (CTEQ3L), PYTHIA 
(CTEQ3L), PYTHIA (CTEQ4L), and PYTHIA (GRV94L) for the integrated 6-quark total cross section {pT> 
5GeV/c, \y\ < 1) at 1.8 TeV. The contributions from flavor creation, flavor excitation, and shower/fragmentation 
are shown together with the overall resulting sum (overall height of box). The DOl.l structure functions include only 
four flavors and the probability of finding a b (or b) quark within the proton is set identically equal to zero. Thus, 
running with the DOl.l structure functions results in no flavor excitation. The differences among the Monte-Carlo 
models for flavor excitation are due to the different ways the models handle the 6-quark mass in this subprocess. 

The Monte-Carlo model predictions for the shower/fragmentation contribution differ considerably. This is not 
surprising since ISAJET uses independent fragmentation, while HERWIG and PYTHIA do not; and HERWIG and 
PYTHIA modify the leading-log picture of parton showers to include "color coherence effects", while ISAJET does 
not. 

The flavor creation contribution to inclusive 6-quark; production {pt> 5GeV/c, \y\ < 1) is predicted by HERWIG 
(CTEQ3L) and PYTHIA (CTEQ3L) to be about 25% of the overaU 6-quark production rate from aU three sources 
at 1.8 TeV, while for ISAJET (CTEQ3L) it is about 35% of the total. 
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FIG. 7: Data on the integrated inclusive &-quark total cross section (pr > pT(min), \y\ < 1) for proton-antiproton collisions at 
1.8 TeV compared with the QCD Monte-Carlo model predictions of HERWIG 5.9 (CTEQ3L, pT(hard) > 3GeV/c). The four 
curves correspond to the contribution from flavor creation (FIG. 0), flavor excitation (FIG. shower/fragmentation (FIG. 
and the resulting total. 
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FIG. 8: The inclusive 6-quark differential cross section, da/dpr, for \y\ <1 for proton-antiproton collisions at 1.8 TeV resulting 
from PYTHIA 6.158 (CTEQ4L, pT(hard)> 0). The four curves correspond to the contribution from flavor creation (FIG. 
flavor excitation (FIG. 0), shower/fragmentation (FIG. 0), and the resulting total. 



IV. CORRELATIONS 

Clearly the three sources of 6-quarks; flavor creation, flavor excitation, and shower/fragmentation have quite different 
topological structures and the correlations between the b and b quarks are quite different. 
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FIG. 9: Data on the integrated inclusive 6-quark total cross section (pr > PT(niin), \y\ < 1) for proton-antiproton collisions at 
1.8 TeV compared with the QCD Monte-Carlo model predictions of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0). The four curves 
correspond to the contribution from flavor creation (FIG. 0), flavor excitation (FIG. shower/fragmentation (FIG. and 
the resulting total. 
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FIG. 10: Predictions of HERWIG 5.9 (DOl.l), HERWIG 5.9 (CTEQ3L), ISAJET 7.32 (CTEQ3L), PYTHIA 6.115 (CTEQ3L), 
PYTHIA 6.158 (CTEQ4L), and PYTHIA 6.115 (GRV94L) for the integrated inclusive 6-quark total cross section (pr > 5 GeV/c, 
\y\ < 1) for proton-antiproton collisions at 1.8 TeV. The contributions from flavor creation (FIG. flavor excitation (FIG. 
shower/fragmentation (FIG. ^ are shown together with the resulting sum (overall height of box). 



A. Simple "Toward" and "Away" Probabilities 



FIGS. |ll| - |13| show the QCD Monte-Carlo model predictions for some simple correlations. Here one requires a 
&-quark to be in the region pT> 5 GeV/c and |y| < 1 and then asks for the probability of finding a &-quark in the same 
region, px> 5 GeV/c and \y \ < 1. Furthermore, one breaks this probability into two terms, "toward" and "away" . The 
"toward" region corresponds to |A(/)|< 90° and the "away" region has |A(/)|> 90°, where A4> is the azimuthal angle 
between the b and b quark. The models predict that for flavor creation the probability that both the b and b quark 
lie in the region pt> 5GeV/c and \y \ < 1 is about 40%, with the 5-quark almost always on the "away-side" of the 
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FIG. 11: Predictions of PYTHIA 6.115 (CTEQ3L, pT(hard) > 3GeV/c) for the probability of finding a &-quark with pr > 
5 GeV/c and \y\ < 1 in an event with a 6-quark with pr > 5 GeV/c and \y\ < 1 for proton-antiproton colhsions at 1.8 TeV. The 
contribution from the "toward" QA(^|< 90°) and the "away" (|A(;/>|> 90°) region of the b-quark are shown for fiavor creation 
(FIG. ^), flavor excitation (FIG. and shower/fragmentation (FIG. ^), where |A(/!>| is the azimuthal angle between the b and 
6- quark. 
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FIG. 12: Predictions of HERWIG 5.9 (CTEQ3L, pT(hard) > 3 GeV/c) for the probability of finding a 6-quark with pr > 5 GeV/c 
and |y I < 1 in an event with a 6-quark with pr > 5 GeV/c and |y | < 1 for proton-antiproton collisions at 1.8 TeV. The contribution 
from the "toward" (|A(^|< 90°) and the "away" (|A(j!)|> 90°) region of the &-quark are shown for flavor creation (FIG. |l|), flavor 
excitation (FIG. H), and shower/fragmentation (FIG. W), where \A(j)\ is the azimuthal angle between the b and fe-quark. 



6-quark. 

The QCD Monte-Carlo models differ considerably on the correlations for the shower/fragmentation contribution. 
For this contribution all three models predict that both the b and b quark lie in the region pT> 5 GeV/c and |y| < 1 
around 30%-40% of the time, which is comparable to the flavor creation contribution. However, ISAJET predicts 
that for the shower/fragmentation contribution that the 6-quark is almost always on the "toward-side" of the 6-quark, 
while HERWIG and PYTHIA predict about equal amounts of "toward" and "away" for this contribution. 

As one would expect, all the models predict that it is not very likely to find both the 6 and the 6-quark in the region 
Pt> 5 GeV/c and |y| < 1 for the flavor excitation contribution. For flavor excitation either the 6-quark or the 6-quark 
is part of the "underlying" event (see FIG. |^). The models differ greatly on the correlations for the flavor excitation 
contribution. For this contribution PYTHIA predicts that both the 6 and 6-quark lie in the region pt> 5 GeV/c and 
|y| < 1 around 20% of the time, while ISAJET predicts 13%_and HERWIG gives 5% for this probability All three 
models predict that for flavor excitation when both the 6 and 6-quark are found in the region pT> 5 GeV/c and |y| < 1 
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FIG. 13: Predictions of ISAJET 7.32 (CTEQ3L, pT(hard) > 3 GeV/c) for the probability of finding a 6-quark with pt > 5 GeV/c 
and < 1 in an event with a 6-quark with pr > 5 GeV/c and |3/| < 1 for proton-antiproton collisions at 1.8 TeV. The contribution 
from the "toward" (|A(;6|< 90°) and the "away" (|A(/)|> 90°) region of the 6-quark are shown for flavor creation (FIG. flavor 
excitation (FIG. ^), and shower/fragmentation (FIG. where \A(I>\ is the azimuthal angle between the b and fe-quark. 

that the 6-quark is about equaUy likely to be on "toward" or "away" side of the 6-quark. 
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FIG. 14: Predictions of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the azimuthal angle, A(j>, between a &-quark with pri > 
5 GeV/c and | j/i | < 1 and a 6-quark with pt2 > and |i/2 1 < 1 in proton-antiproton collisions at 1.8 TeV. The curves correspond 
to da/dAfj) (/ib/°) for flavor creation (FIG. ^), flavor excitation (FIG. ^), shower/fragmentation (FIG. ^, and the resulting 
total. The "toward" (|A(;/)|< 90°) and the "away" (|A(j!)|> 90°) region of the 6-quark are labeled. {Note the logarithmetic scale.) 



B. A0 Distribution 

FIG. 1^ shows the predictions of PYTHIA for the azimuthal angle, A0, between a 6-quark with pn > 5 GeV/c 
and \yi\ < 1 and b quark with pT2 > and I2/2I < 1- FIG. |l^ is the same as FIG. |lj except both the 6 and 6 quark 
are required to have px > 5GeV/c. Notice that the relative amounts of flavor excitation and shower/fragmentation 
decrease compared to flavor creation when one demands that both the 6 and 6 quarks to have pt > 5GeV/c. 
Integrating the curves in FIG. ^ over the "toward" and "away" regions correspond to the simple probabilities 
discussed in Section IV(a). Clearly, the "toward" region is very sensitive to the presence of the flavor excitation and 
shower/fragmentation terms. 
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FIG. 15: Predictions of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the azimuthal angle, A(j>, between a &-quark with pri > 
5GeV/c and \yi\ < 1 and a fe-quark with pt2 > 5GeV/c and \y2\ < 1 in proton-antinroton collisions at 1.8 TeV. The curves 
correspond to da/dA(f> (/ib/°) for flavor creation (FIG. flavor excitation (FIG. shower/fragmentation (FIG. and 
the resulting total. The "toward" (|A(^|< 90°) and the "away" (|A</)|> 90°) region of the 6-quark are labeled. {Note the 
logarithmetic scale.) 
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FIG. 16: Predictions of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the distance, _R, in rj-(\> space between a &-quark with 
Pti > 5GeV/c and \y\\ < 1 and a 6-quark with pt2 > 5GeV/c proton-antiproton collisions at 1.8 TeV. The curves correspond 
to da/dR (/ib) for flavor creation (FIG. 0), flavor excitation (FIG. shower/fragmentation (FIG. and the resulting total. 
(Note the logarithmetic scale.) 



C. Distance "R" 

FIG. ^ shows the predictions of PYTHIA for the distance, R, in rj-cj) space between a 6-quark with pTi > 5 GeV/c 
and \yi\ < 1 and 6-quark with pT2 > 5GeV/c, where R = ^ (r\\ — 772)^ -I- (0i — 02)^ and 77 is the pseudo-rapidity. 
One can see that the region R < \ isolates the shower/fragmentation contribution. 

D. Rapidity Correlations 

FIG. |l^ shows the predictions of PYTHIA for the rapidity, yi of a 6-quark with p^i > 5 GeV/c for events with a 
6-quark with 'pT2 > 5 GeV/c and \y2\ < 0.5. Here the flavor excitation contribution behaves much differently than the 
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other two terms. However, the resulting sum of the three contributions looks similar to the flavor creation term. 
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FIG. 17: Predictions of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the rapidity, yi, of a 6-quark with pri > 5GeV/c for 
events with a 6-quark with pt2 > 5GeV/c and I2/2I < 0.5 in proton-antiproton colhsions at 1.8 TeV. The curves correspond to 
da/dyi (fih) for flavor creation (FIG. |lj), flavor excitation (FIG. ^), shower/fragmentation (FIG. ^), and the resulting total. 
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FIG. 18: Prediction of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the pt asymmetry, A = (pri -Pt2)/{pti +PT2), for events 
with a &-quark with pri > and < 1.0 and a 6-quark with pt2 > 5 GeV/c and \y2\ < 1.0 in proton-antiproton collisions at 
1.8 TeV. The curves correspond to da/dA {ph) for flavor creation (FIG. 0), flavor excitation (FIG. ^), shower/fragmentation 
(FIG. 0), and the resulting total. 



E. The Transverse Momentum Asymmetry 

FIG. shows the predictions of PYTHIA for the px asymmetry, A — {pxi — Pt2)/{pti +PT2), for events with 
a 6-quark with p^i > and |yi| < 1 and a 6-quark with pt2 > 5GeV/c and \y2\ < 1. If there were no initial- 
state radiation (see FIG. |^) and no "intrinsic" transverse momentum, the flavor creation terms would produce zero 
asymmetry. Also, as one would expect (FIG. the flavor excitation terms have, on the average, a large transverse 
momentum asymmetry. Requiring px asymmetries with a magnitude greater than around 60% isolates the flavor 
excitation contribution. 
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FIG. 19: Prediction of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the transverse momentum, pt2, of a 6-quark with \y2\ < 1.0 
for events with a b-quark with pri > 5 GeV/c and \yi \ < 1.0 in proton-antiproton coUisions at 1.8 TeV. The curves correspond to 
da/dpT2 (/ib/GeV/c) for flavor creation (FIG. flavor excitation (FIG.^), shower/fragmentation (FIG. and the resulting 
total. {Note the logarithmetic scale.) 
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FIG. 20: Prediction of PYTHIA 6.158 (CTEQ4L, pT(hard)> 0) for the integrated pair cross section for a 6-quark with 
PT2 > pT2(min) and \y2\ < 1.0 and a 6-quark with pri > 5GeV/c and \yi\ < 1.0 in nroton-antiproton collisions at l.STeV. 
The curves correspond to a (/ib) for flavor creation (FIG. flavor excitation (FIG. shower/fragmentation (FIG. and 
the resulting total. {Note the logarithmetic scale.) 



F. Double Differential Cross Section 



FIG. ^ shows the predictions of PYTHIA for the transverse momentum, pt2, of a &-quark with \y2\ < 1 for events 
with a b-quark with pxi > 5GeV/c and \yi\ < 1. As the transverse momentum asymmetry in FIG. ^ indicated, 
the flavor creation contribution peaks at pt2 ~ 5GeV/c {i.e., small pr asymmetry), while the flavor excitation and 
shower/fragmentation contributions peak at low (i.e., larger asymmetry). 

FIG. shows the predictions of PYTHIA for the integrated bb pair cross section for a &-quark with pT2 > PT2(inin) 



13 



and \y2\ < 1 and b-quark with pxi > 5GeV/c and \yi\ < 1. Notice how the three terms contribute in different 
proportions to the bb pair cross section in FIG. ^ and the inclusive single 5-quark cross section in FIG. ||. Flavor 
excitation is a major contributor to the single inclusive cross section, while flavor creation dominates the pair cross 
section. 



V. SUMMARY AND CONCLUSIONS 

The leading-log order QCD hard scattering Monte-Carlo models of HERWIG, ISAJET, and PYTHIA have been 
used to study the sources of &-quarks at the Tevatron. The reactions responsible of producing fo-quarks are separated 
into flavor creation, flavor excitation, and shower/fragmentation. 

Flavor Creation 

The production of a bb pair via gluon fusion or annihilation of light quarks is easy to generate and all three 
QCD Monte-Carlo models {with the same structure functions) predict roughly the same cross section and similar 
bb correlations from these terms. At the Tevatron all three Monte-Carlo models predict that the flavor creation 
contribution to 6-quark production is less than 35% of the overall inclusive b-quark production rate {pt> 5GeV/c, 
|t/| < 1) from all three sources. 

Flavor Excitation 

The source of &-quarks resulting from the scattering of a 6-quark or &-quark out of the initial-state into the final-state 
by a gluon or by a light quark or antiquark is difficult to generate and depends sensitively on the parton distribution 
functions. The QCD Monte-Carlo model predictions differ somewhat, however, it seems likely that at the Tevatron the 
flavor excitation contribution to the inclusive 6-quark cross section is comparable to or greater than the contribution 
from flavor creation. The bb correlations resulting from the flavor excitation term are much different than those arising 
from flavor creation. In particular, the distribution in the azimuthal angle, A0, between the b and the 6-quark peaks 
sharply in the "backward" (|A0|> 90°) region for flavor creation, while flavor excitation produces an approximately 
flat |A(/)| distribution. Also, the flavor excitation terms have, on the average, a large transverse momentum asymmetry 
between the b and &-quark, while the flavor creation terms produce a pT asymmetry near zero. 

Shower/Fragmentation 

The production of bb pairs within a parton-shower or during the fragmentation process of gluons and light quarks 
and antiquarks is an important source at the Tevatron. The QCD Monte-Carlo models predictions differ considerably 
for this contribution. However, at the Tevatron the shower/fragmentation contribution to the inclusive b-quark cross 
section might be comparable to the contribution from flavor creation. This source can be isolated by looking for bb 
pairs with i? < 1, where R is the distance in r/-0 space between the 6 and 6-quark. 

The QCD Monte-Carlo model leading-log estimates of the flavor excitation and the shower/fragmentation contri- 
butions to 6-quark production are uncertain and should not be taken too seriously. However, it seems likely that at 
the Tevatron all three sources of 6-quarks, flavor creation, flavor excitation, and shower/fragmentation are important 
One does not expect precise agreement from leading-log estimates. On the other hand, the qualitative agreement 
shown indicates that probably nothing unusual is happening in 6-quark production at the Tevatron. Furthermore, in 
Run II at the Tevatron we should be able to isolate experimentally the individual contributions to 6-quark production 
by studying 66 correlations in detail. 
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